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ABSTRACT 

Using a redshift survey of 1323 galaxies (1092 new or remeasured) in a region 
of 95 square degrees centered on the nearby galaxy cluster Abell 2199, we ana- 
lyze the supercluster containing A2199, A2197, and an X-ray group. The caustic 
technique accurately reproduces the true mass profiles of simulated simple super- 
clusters (i.e., superclusters where the virial mass of one cluster is 2-10 times the 
virial mass of all other clusters in the supercluster). We calculate the masses of 
the two main components of A2197 (A2197W and A2197E) using archival X-ray 
observations and demonstrate that the A2199 supercluster is simple (the mass of 
A2199 is 5 and 12 times larger than A2197W and A2197E respectively) and thus 
that the caustic technique should yield an accurate mass profile. The masses of 
A2199, A2197W, and A2197E within rsoo (the radius within which the enclosed 
density is 500 times the critical density) are 22.0, 3.8, and 1.7xlO 13 /i _1 M . The 
mass profile is uncertain by ~30% within 3 /i _1 Mpc and by a factor of two within 
8 /i _1 Mpc and is one of only a few for a supercluster on such large scales. In- 
dependent X-ray mass estimates agree with our results at all radii where they 
overlap. The mass profile strongly disagrees with an isothermal sphere profile 
but agrees with profiles suggested by simulations. We discuss the interplay of 
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the supercluster dynamics and the dynamics of the bound subclusters. The agree- 
ment between the infall mass profile and other techniques shows that the caustic 
technique is surprisingly robust for simple superclusters. 

Subject headings: galaxies: clusters: individual (A2199, A2197) — galaxies: kine- 
matics and dynamics — cosmology: observations 

1. Introduction 

Clusters of galaxies are important probes of the distribution of matter on intermediate 
scales. Clusters are surrounded by infall regions where the galaxies are falling into the 
gravitational potential well of the cluster, but have not yet reached equilibrium. Many, 
perhaps most, of the galaxies in this region are on their first pass through the cluster. They 
populate a regime between that of relaxed cluster cores and the surrounding large-scale 
structure where the transition from linear to non-linear clustering occurs. Recent papers 
have explored infall regions using two-body dynamics of binary clusters (e.g., Mohr & 
Wegner 1997), the virial theorem in superclusters (Small et al. 1998), and weak lensing 
(Kaiser et al. 2002). Ellingson et al. (2001) demonstrate that these infall regions are the site 
of significant galaxy evolution and are key to understanding the Butcher-Oemler (Butcher 
& Oemler 1984) effect. Two recent papers (Reisenegger et al. 2000; Rines et al. 2001b) 
demonstrate that the infall regions of clusters (the Shapley Supercluster and Abell 2199 
respectively) provide evidence of hierarchical structure formation on supercluster scales. 

Because clusters are not in equilibrium outside the virial radius, neither X-ray observa- 
tions nor virial analysis provide accurate mass determinations at large radii. There are now 
two methods of approaching this problem: weak gravitational lensing (Kaiser et al. 2002) 
and kinematics of the infall region (Diaferio & Geller 1997; Diaferio 1999, hereafter DG 
and D99 respectively). In redshift space, the infall regions of clusters form a characteristic 
trumpet-shaped pattern. These caustics arise because galaxies fall into the cluster as the 
cluster potential overwhelms the Hubble flow (Kaiser 1987; Regos & Geller 1989). Under 
simple spherical infall, the galaxy phase space density becomes infinite at the caustics. DG 
analyzed the dynamics of infall regions with numerical simulations and found that in the 
outskirts of clusters, random motions due to substructure and non-radial motions make a 
substantial contribution to the amplitude of the caustics which delineate the infall regions 
(see also Vedel & Hartwick 1998). DG showed that the amplitude of the caustics is a mea- 
sure of the escape velocity from the cluster; identification of the caustics therefore allows 
a determination of the mass profile of the cluster beyond the virial radius. DG find that 
nonparametric measurements of caustics in simulated clusters yield mass profiles accurate 
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to ~50% on scales up to 10 h 1 Mpc. 

The caustic technique places no requirements on the equilibrium state of the cluster; 
rather, it assumes only that galaxies trace the velocity field. Indeed, simulations suggest 
that little or no velocity bias exists on linear and mildly non-linear scales (Kauffmann et al. 
1999a,b). Vedel & Hartwick (1998) used simulations to explore an alternative parametric 
analysis of the infall region using a maximum likelihood method. Their technique requires 
assumptions about the functional forms of the density profile and the velocity dispersion 
profile. Geller et al. (1999, hereafter GDK) applied the kinematic method of DG to the 
infall region of the Coma cluster. They successfully reproduced the X-ray derived mass 
profile and extended direct determinations of the mass profile to ~10 h~ 1 Mpc. This method 
has also been applied to the Shapley Supercluster (Reisenegger et al. 2000), A576 (Rines 
et al. 2000), A1644 (Tustin et al. 2001) and the Fornax cluster (Drinkwater et al. 2001). 

We apply the method of DG and D99 to A2199, one of the richest, most regular nearby 
clusters of galaxies (Markevitch et al. 1999, hereafter MVFS). Rines et al. (2001b) presented 
a preliminary view of this system and demonstrated that it contains several infalling, bound 
subclusters. This system provides an important test of the robustness of the caustic technique 
in simple superclusters (i.e., those dominated by a cluster 2-10 times more massive than any 
other system in the supercluster). 

We describe numerical simulations of simple superclusters in § 2 to demonstrate the 
validity of the caustic technique for these systems. We describe the spectroscopic observa- 
tions in § 3. In § 4, we show that A2197 is a superposition of groups and demonstrate that 
A2199/A2197 is indeed a simple supercluster. We determine the amplitude of the caustics, 
calculate the mass profile, and compare the infall mass profile to other estimators in §5. We 
discuss our results in §6 and conclude in §7. The physical scale at the redshift of the super- 
cluster (cz = 9156 km s" 1 ) is 1° = 1.54 h^Mpc (H = 100 h km s~\ ft m = 0.3, fi A = 0.7). 

2. Simulations of Simple Superclusters 

The caustic technique as developed by D99 is successful in reproducing the mass profiles 
of isolated clusters. Before applying the technique to more complex systems, we need to 
verify that this technique still produces accurate mass profiles in the presence of significant 
substructure in the infall region. We extracted three simple superclusters from the ACDM 
cosmological simulation of the GIF collaboration, where galaxy formation and evolution 
within dark matter halos are modelled with a semi-analytic technique (Kauffmann et al. 
1999a,b). Here, we define simple superclusters as those dominated by a cluster 2-10 times 
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more massive than any other systems in the supercluster. These three systems have mass, 
Mi, within r 2 oo (roughly equivalent to the virial mass) for the central system greater than 
1O 14 /i _1 M ; the next largest subhalo has a mass M 2 between 1O 13 /i" 1 M and 1O 14 /i _1 M and 
lies outside the virial radius of the main cluster. Table 1 lists the properties of the simulated 
superclusters. Column [1] gives the system ID, Column [2] is the mass Mi of the largest 
system in the supercluster, Column [3] is the mass ratio Mi/M 2 , Column [4] and Column 
[5] list the mass and 3-D distance of the four subhalos (with M > 10 13 /i _1 M o ) closest to the 
largest system. The mass ratio M x jM 2 ranges between 2 and 6 to investigate the sensitivity 
of the mass estimator to companions of various relative sizes. We restrict our analysis to 
ACDM simulations because these more accurately reproduce the appearance of the observed 
caustics in nearby systems (see references in §1); caustics in rCDM simulations are less 
well-defined. D99 suggests that the degree of delineation of the caustics is a cosmological 
indicator; caustics are better defined (have fewer interlopers) in low-Q m universes because 
the accretion rate in the present epoch decreases with decreasing Q m . 

Figures 1-3 display the caustics (calculated according to the prescription of D99) for 
these systems as viewed from three different lines of sight. Redshift diagrams include all 
galaxies brighter than Mb = —18.5 + 5\ogh. As in isolated clusters (D99), some lines of 
sight are more favorable than others. Figure 4 compares the mass profiles inferred from the 
caustics to the true radial mass profiles. Figure 5 shows the difference between the estimated 
and true mass profiles in units of the uncertainty in the estimated mass. The estimated mass 
differs by more than 2-er only for some lines of sight to the supercluster with M 1 /M 2 = 2.0. 
The bias in the mass estimates appears to be small, although the small number of systems 
considered prevents firm conclusions. Thus, the caustic method is robust to the presence of 
massive subhalos provided the mass ratio (measured within r2oo) is greater than about 2.0. 
We show in §4 that A2197 (actually A2197W+A2197E) is much less massive than A2199, 
the caustic technique should thus yield an accurate mass profile for the A2199 supercluster. 



2.1. Model Mass Profiles 

We now investigate whether the mass profiles of the simulated simple superclusters are 
well described by simple analytic models. The simplest model of a self-gravitating system 
is a singular isothermal sphere (SIS). The mass of the SIS increases linearly with radius. 
Navarro et al. (1997) and Hernquist (1990) propose two-parameter models based on CDM 
simulations of haloes. At large radii, the NFW mass profile increases as ln(r) and the mass 
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of the Hernquist model converges. The NFW mass profile is 

where a is the scale radius and M(a) is the mass within a. We fit the parameter M(a) 
rather than the characteristic density 8 C (M(a) = 47r<5 c p c a 3 [ln(2) — |] where p c is the critical 
density) because M(a) and a are much less correlated than 5 C and a. The Hernquist mass 
profile is 

M « '') = M (7T^ (2) 

where a is the scale radius and M is the total mass. Note that M(a) = M/A. The SIS mass 
profile is 

M« r) = M(a = 0.5 fc -. M pc) jjJ ^ (3) 

where we arbitrarily set the scale radius a = 0.5 h^Mpc. We assume 10% uncertainties in 
the true mass profiles and minimize x 2 - Table 2 lists the best-fit parameters a (fixed for 
SIS) and M(a) for the three models for the simulated superclusters. Because the individual 
points in the mass profile are not independent, the absolute values of xl listed in Table 2 are 
indicative only, but it is clear that the NFW and Hernquist profiles provide acceptable fits 
to the infall mass profile; the SIS is excluded for all estimates. The NFW profile generally 
provides a better fit to the data than the Hernquist profile. A non-singular isothermal sphere 
mass profile yields results similar to the SIS; thus, we only report our results for the SIS. 



3. Observations 
3.1. Spectroscopy 

We have collected 1323 redshifts (1092 new or remeasured) in a large region (95 square 
degrees) surrounding A2199/A2197. Two of us (JJM and GW) used the Decaspec (Fabricant 
& Hertz 1990) at the 2.4-m MDM telescope on Kitt Peak to obtain 246 redshifts of galaxies 
in the central 1° x 1° of this region for a separate Jeans analysis of the central region of 
A2199 (Mohr et al. in preparation). 

We used the FAST spectrograph (Fabricant et al. 1998) on the 1.5-m Tillinghast tele- 
scope of the Fred Lawrence Whipple Observatory (FLWO) to measure 820 spectra of galaxies 
within 5°.5 («8.5 h _1 Mpc) of the center of A2199. FAST is a high throughput, long slit 
spectrograph with a thinned, backside illuminated, antireflection coated CCD detector. The 
slit length is 180"; our observations used a slit width of 3" and a 300 lines mm" 1 grating. 
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This setup yields spectral resolution of 6-8 A and covers the wavelength range 3600-7200 
A. We obtain redshifts by cross-correlation with spectral templates of emission-dominated 
and absorption-dominated galaxy spectra created from FAST observations (Kurtz & Mink 
1998). 

We observed infall galaxy candidates in three campaigns. We selected targets from dig- 
itized images of the POSS I. We initially selected galaxies using the automatic classification 
system of the Automated Plate Scanner (APS) 8 ; we visually inspected these targets to elim- 
inate stars. The first campaign yielded a deep sample in the central 4° x 4° region around 
A2199. This sample is complete to 103aE magnitude E < 16.5 and consists of 304 redshifts. 
The second campaign (379 redshifts) is a shallower survey (E < 16.1) of all galaxies within 
5.°5 ph 8.5 /i -1 Mpc of A2199. The third campaign added 137 redshifts for galaxies which 
were excluded from the APS catalog for unknown reasons (§3.2). The completeness limits 
are imprecise because the magnitudes come from multiple plate scans and because we could 
not obtain redshifts for some low surface brightness galaxies. We include 84 redshifts asso- 
ciated with the groups NRGs385 and NRGs388 obtained with FAST for a separate study of 
the X-ray and optical properties of groups of galaxies (59 published in Mahdavi et al. 1999, 
2000, hereafter MBGR). 

We collected the remaining redshifts from the Updated Zwicky Catalogue (Falco et al. 
1999, 158 redshifts) with a small number (15) from other sources listed in the NASA/IPAC 
Extragalactic Database 9 (Kirshner et al. 1983; Strauss & Huchra 1988; de Vaucouleurs et al. 
1991; Freudling et al. 1992; Zabludoff et al. 1993; Haynes et al. 1997; Hill & Oegerle 1998). 
Table 3 contains our redshift data. The table includes all redshifts obtained with FAST 
and from the literature as well as redshifts obtained with Decaspec for all galaxies brighter 
than £=16.1 (the remaining 188 will appear in Mohr et al. in preparation). This catalog is 
complete to £=16.1 within 5°. 5 of the center of A2199 and includes 38 galaxies outside 5°. 5 
but within 6°. 5 w 10 /i _1 Mpc (Figure 6, described in detail in §5). 

3.2. APS Catalog - Photometry and Incompleteness 

When matching galaxies from our redshift catalog with a galaxy catalog from APS, we 
discovered that several nearby, bright galaxies including NGC 6160 (a cD galaxy in A2197) 
are not classified as galaxies by the APS neural network. A further search reveals that the 

8 The APS databases are supported by the National Aeronautics and Space Administration and the 
University of Minnesota, and are available at http://aps.umn.edu 

9 The NASA/IPAC Extragalactic Database is available at http://nedwww.ipac.caltech.edu/index.html 
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APS identified many of these objects but incorrectly classified them as stars. The default 
magnitude listed for these objects is from the magnitude-diameter relationship for stars 
(Humphreys et al. 1991). Because galaxies have large diameters, they are assigned bright 
magnitudes. We extracted a catalog of all objects classified as stars by APS with magnitudes 
Ediam < 12 within 5.°5 of A2199. Visual inspection of these 4876 objects yielded 160 galaxies 
(5 of which are star-galaxy blends), suggesting a misclassification rate of ~3.3%. Because 
these objects include bright, large galaxies like NGC 6160, the incompleteness in the APS 
could be significant in nearby galaxy surveys. 

Because of these difficulties, we reanalyzed the POSS 1 103aE (red) plates using SExtrac- 
tor to calculate object parameters and magnitudes. We calculated the mean offset between 
SExtractor MAGBEST magnitudes and the isophotal magnitudes in the APS catalog to de- 
termine the zero point for each plate. We next identified galaxy candidates as objects with 
MAGBEST<16.5, staricity parameter smaller than 0.5, and IS07/IS02 greater than 0.35. 
ISO 1-7 are the numbers of pixels included within various isophotes; the latter criterion thus 
selects extended objects. We ran SExtractor in CCD mode. We experimented with running 
SExtractor in PHOTO mode (appropriate for photographic plates) after taking the loga- 
rithm of the scanned image (this step is needed to obtain a sufficiently small dynamic range 
for this mode) and found nearly identical results for both photometric and shape parameters. 

Finally, we visually inspected all galaxy candidates to compile a complete catalog of 
galaxies brighter than E — 16.1. Our new catalog includes 970 galaxies, 135 (14%) of which 
are not included in the APS catalog. These 970 galaxies comprise the majority (71%) of the 
redshift sample; the fainter galaxies in the remainder of the sample are located mainly near 
A2199, A2197, and the X-ray groups NRGs385 and NRGs388. 

3.3. X-Ray Data 

Archival X-ray observations are available for the entire A2199 supercluster from the 
ROSAT All-Sky Survey (Voges et al. 1999). Several longer pointed observations of selected 
parts of the supercluster are available in the ROSAT PSPC and ASCA archives (Table 4). 
Feretti et al. (1995) discuss radio and ROSAT pointed observations of the groups NRGs385 
and NRGs388. We calculate the luminosities of the systems in the ROSAT energy band (Ta- 
ble 5) from the pointed observations assuming that all of the subsystems are at the distance 
of A2199 (cz =9156 km s"\ see §5). We used the X-ray analysis package ZHTOOLS 10 to 
extract regions around the sources, subtract background sources, and correct for vignetting 



10 Available at http://hea-www.harvard.edu/ alexey/zhtools/ 
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and the exposure maps. We then used PIMMS 11 to convert PSPC count rates into fluxes. 
We assumed Galactic column density and a Raymond-Smith thermal plasma with an abun- 
dance of 0.4 and kT= 4.5 keV for A2199 (MVFS), 0.2/1.7 keV for A2197W, and 0.2/1.0 
keV for all others (see §5.2). Figure 6 displays the RASS data in a large area around A2199 
after correcting for the exposure map and smoothing with a Gaussian with a FWHM of 2'. 

4. A2197: An Optically Rich "Cluster" Composed of Superposed Groups 

A2197 has been observed by both RO SAT and ASCA. Figure 7 shows X-ray contours of 
the 13.5-ksec combined archival RO SAT PSPC observation of A2197 (Table 4) overlaid on 
an optical image from POSS I. The X-ray contours clearly show the presence of West-East 
components separated by 24' centered on NGC 6160 and NGC 6173 respectively (Muriel 
et al. 1996). The galaxy contours (Figure 6) show at least two distinct NW-SE components, 
but the galaxy distribution is skewed with respect to the X-ray emission (see also Rines et al. 
2001b). Because the intragroup gas is expected to better trace the gravitational potential, 
we use the X-ray centers in all subsequent analysis. 

We find several distinct X-ray sources in the combined PSPC observation of A2197. 
Table 6 lists the X-ray luminosities of sources with obvious optical counterparts. In addition 
to A2197W and A2197E, we identify five galaxies in or near A2197, one background galaxy, 
two background Abell clusters (A2187 and A2196), and two background QSOs. A2187 is 
recognized as an X-ray cluster (Ebeling et al. 2000; Bohringer et al. 2000) but our deter- 
mination of the X-ray luminosity of A2196 is new. The galaxies have X-ray luminosities of 
L x (0A — 2.4 keV) = 4 — 64 x 10 40 /i~ 2 erg s -1 , within the range of isolated bright galaxies 
(Mahdavi & Geller 2001). 

4.1. Redshift Distribution 

A2197 shows significant substructure in redshift space. Figure 8 shows a redshift his- 
togram of galaxies within 1.5 /i _1 Mpc of A2197W with a bin size of 200 km s" 1 . We cal- 
culate x 2 = 36.3 for 12 degrees of freedom, indicating that the parent redshift distribu- 
tion is non-Gaussian at the 99.9% confidence level. A two-component Gaussian fit yields 
X 2 = 10.8 for 9 degrees of freedom, an acceptable fit. The properties of the best-fit compo- 
nents are Ni = 82 ± 5, ^ = 8787 ± 33 km s~\ oi = 361 ± 26 km s~\ and N 2 = 48 ± 6, 



Available at http://heasarc.gsfc.nasa.gov/docs/software/tools/pimms.html 
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/i 2 = 9800 ± 35 km s"\ a 2 = 365 ± 78 km s -1 , 

We recalculate the mean redshifts and velocity dispersions of the two subclusters ac- 
cording to the prescription of Danese et al. (1980) but with a limiting radius of 0.25 /i -1 Mpc, 
which is small enough that the groups do not overlap. The mean redshifts of A2197W and 
A2197E are 9569 ± 151 km s _1 and 8919 ± 153 km s' 1 based on 11 and 10 members respec- 
tively. The mean redshifts differ by 4.2-cr. A K-S test indicates that the redshifts are drawn 
from different populations at the 99.6% confidence level. Note that the redshifts of NGC 
6160 and NGC 6173 are 9408 km s _1 and 8842 km s _1 respectively (Table 3), suggesting 
that they may lie near the centers of the two groups (Figure 8). 

These tests show that the groups A2197W and A2197E are separated in redshift space. 
The two groups overlap each other both on the sky and in redshift space; they may be 
interacting in the current epoch. Apparently, A2197E is a richer group than A2197W but 
has a lower X-ray luminosity. The superposition of the two groups artificially enhances the 
apparent optical richness of the system. 

4.2. Surface Brightness Distribution 

We fit a f3 model (Cavaliere & Fusco-Femiano 1976) convolved with the PSPC PSF to 
the observed surface brightness of A2197W in the PSPC observation. We find a core radius 
r c = 0J21 = 11 h^kpc and (3 X = 0.41. Note that A2197E is 28' off-axis, the PSF of the 
PSPC at this radius is ~1' and is anisotropic. Because of these issues, we do not attempt to 
fit the surface brightness profile of A2197E. Our results agree with Muriel et al. (1996). 

4.3. ICM Temperature 

ASCA's broad energy band (0.5-10.0 keV) is particularly useful for determining cluster 
temperatures. Because of the poor angular resolution of ASCA, we determine the emission 
weighted average temperature within 11', or ~ 0.3 h^Mpc; outside this radius the cluster 
emission is overwhelmed by the background. We obtained the screened data from Goddard 
Space Flight Center (GSFC). We extract a spectrum including all photons within a circle of 
radius 11' (44 pixels) centered on the cluster center for GIS data from long observations of 
A2197W and A2197E (Table 4). The centroid of the SIS images agree with the positions of 
the ROSAT centroids within the ~ 0^4 uncertainty in the SIS position (Gotthelf 1996) 12 . 



Available at http://heasarc.gsfc.nasa.gov/docs/asca/newsletters/Contents4.html 
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Using XSPEC (vlO.O), we fit the cluster spectrum to a model including absorption 
('wabs') parameterized by the column density of hydrogen (which we set to the Galactic 
value) and the standard Raymond-Smith model (Raymond & Smith 1977) characterized by 
temperature, iron abundance, redshift, and a normalization factor. The iron abundance 
is measured relative to cosmic abundance. We fit the temperature, iron abundance, and 
normalization as free parameters. Because we are interested in the emission-weighted X-ray 
temperature, we use only the GIS data because the SIS field of view does not contain the 
entire X-ray emitting regions. 

We use the weighting system developed by Churazov et al. (1996), which properly 
accounts for errors in small-number statistics. Because there are few counts above 8.0 keV 
and the calibration uncertainties become large at low energies, we include only data from 
0.8 — > 8.0 keV. We fit the spectrum to slightly different ranges to ensure that the fitted 
model parameters are consistent for different choices. 

We obtain acceptable fits by assuming that the gas is isothermal. More complicated 
models are thus unnecessary. Our best-fit model for A2197W has an ICM temperature of 
1.55+°io keV with an iron abundance 0.13t°;^ cosmic. For A2197E, we find kT '=0.961°;^ 
keV and an iron abundance 0.20lg;^ cosmic. Uncertainties are 90% confidence limits for one 
parameter and do not include calibration uncertainties. The acceptable isothermal model 
fits are consistent with no large temperature variations within either A2197W or A2197E, 
suggesting an absence of shock fronts which should be present if the groups were in the 
process of merging. Figure 9 (10) shows the X-ray spectrum and the best-fit model for 
A2197W (A2197E). These temperatures are slightly smaller but consistent with those in 
Finoguenov et al. (2001a) estimated from ASCA SIS data for the central regions of the 
clusters. 

Using the best-fit model, we calculate a bolometric flux of 5.1(4.3) x 10~ 12 ergs cm~ 2 
s" 1 for A2197W (A2197E). We calculate a bolometric luminosity of 5.3(4.5) x 10 42 h~ 2 ergs 
s _1 . If we use the ASCA temperatures and abundances and the ROSAT PSPC count rates 
from the pointed observations, we estimate bolometric luminosities of 5.0(4.3) x 10 42 /i~ 2 ergs 
s -1 , in good agreement with the ASCA estimates. 

4.4. Mass Estimates for A2197W and A2197E 

The mass-temperature relation (Horner et al. 1999; Nevalainen et al. 2000; Finoguenov 
et al. 2001b) gives a straightforward estimate of the mass of a cluster given its temperature. 
Finoguenov et al. (2001b) find a relation of M 500 = (1.87 ± 0.14) x lO^T^^h^M^ 
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where M 50 o is the enclosed mass within a radius r^o inside which the mean density is 500 
times the critical density and T^y is the electron temperature in keV. Using T} ce y=4.5 
(MVFS) for A2199, the values of M 500 for A2199, A2197W, and A2197E are 22.0, 3.8, and 
1.7 respectively (in units of 10 13 /i _1 M Q ). Thus, the mass of A2199 is ~4 times larger than 
the combined masses of A2197W and A2197E. The A2199 system is therefore similar to the 
systems considered in §2 so we are justified in applying the caustic technique to A2199. 

With the assumptions of hydrostatic equilibrium, negligible non-thermal pressure, and 
spherical symmetry, the gravitational mass inside a radius r is 

, , , n kT /d\np aas dlnT\ . A , 

M « (<r) = -^hnff + dl^) r (4 > 

(Fabricant et al. 1980). For a uniform temperature distribution, the second term on the right 
hand side vanishes. We then only need to determine the gas temperature and the density 
distribution of the gas to calculate the gravitational mass. Under the standard hydrostatic- 
isothermal j3 x model, the mass is related to j3 x and a by 

M « (< r) = l^rhm = 5 ' 65 X IG^eV^-'Me (5) 

where M to t{< f) is the total gravitational mass within a radius r and the numerical approx- 
imation is valid for T^ e y in keV and r and a in /i _1 Mpc. 

We use Equation 4 to calculate the total mass of A2197W within 0.5 h~ 1 Mpc; the 
estimate outside ~0.3 /i _1 Mpc is an extrapolation. We find a mass of M ~ 2.2 and 3.6 x 
10 13 /i _1 M at r = 0.3 and 0.5 h~ 1 Mpc respectively. White et al. (1997) use spherical 
deprojection to estimate M ~ 3.9 x 10 13 h~ l M & at r = 0.28 /i -1 Mpc for A2197W, in 
reasonable agreement with our results. 

We calculate the virial mass of A2197 from galaxies within 1.5 /i _1 Mpc of the X-ray 
peak of A2197W. We find M vir = 5.4 x 10 14 /i" 1 M Q , significantly larger than the X -ray mass 
extrapolated assuming an isothermal sphere. Girardi et al. (1998) similarly find that the 
virial mass of A2197 is larger than the X-ray mass. This discrepancy is easily understood 
with the recognition that A2197 is the superposition of at least two groups. Supercluster 
galaxies not contained in the groups may further contribute to an artificially large observed 
velocity dispersion and hence virial mass. The calculated virial mass for A2197 is therefore 
meaningless. 
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5. Defining the Infall Region with Caustics 

Figure 11 displays the redshifts of galaxies surrounding A2199 as a function of projected 
radius. The A2199 supercluster is located within the Great Wall (Geller & Huchra 1989, 
Figure 6b of Falco et al. 1999), which may complicate the interpretation of the dynamics of 
the system. The expected caustic pattern is easily visible. The caustic diagram of A2199 
shows fewer interlopers than in the simulated superclusters, perhaps due to the galaxy for- 
mation recipe used in the simulations or perhaps due to lower infall rates such as would be 
present in a cosmology with Vt m < 0.3. 

We perform a hierarchical structure analysis to locate the centroid of the largest system 
in the supercluster (D99). This analysis yields a position for A2199 of a = 16^28 m 47!0, 5 = 
39°30'22" (J2000) and cz = 9156 km s -1 . This position is 3fl = 80 /i^kpc SE of the X-ray 
peak. 

Rines et al. (2001b) showed that the sky positions of galaxies in the infall region reveal 
several groups in addition to the main cluster. Figure 6 displays X-ray intensity from the 
RASS with contours of the local galaxy density overlaid. We calculate the galaxy density 
contours from cluster members (see definition in §5.7). Extended X-ray emission usually 
comes from an intragroup medium and thus confirms that many, but not all, of the systems 
are bound. NRGs395, NRGs399, and NRGs400 are X-ray faint (L x < 10 42 /i- 2 ergs s" 1 ) 
and thus may not be bound systems; we drop these from further analysis. Table 5 lists 
the coordinates and basic properties of the bound systems (see also MBGR). We calculate 
redshift centers cz and projected velocity dispersions a p according to Danese et al. (1980) 
for galaxies within 1.5 h^Mpc of the X-ray positions. Table 5 lists these properties and the 
number N 3ap of member galaxies (within ±3cx of the central redshift) from which they are 
calculated. Note that the mean redshift of A2199 is consistent with the hierarchical center, 
which we use in all other analyses. A2197W and A2197E overlap both on the sky and in 
redshift space (see §4). We caution the reader that the calculated values of a p for the groups 
may be unreliable estimates of their dynamical properties (see §6.1). 

The region around NGC 6159 shows a density enhancement in both X-rays and galaxy 
density (visible in Figure 2 of Rines et al. 2001b, but not discussed there); A2192, a back- 
ground cluster, lies immediately to the West. Trinchieri & Pietsch (2000) discuss a deep HRI 
observation of NGC 6159 (targeted for its unusually high ratio of X-ray to optical flux); the 
emission is extended. They suggest that the emission is intragroup emission but they lack 
the redshift information necessary to confirm this suggestion. Our redshift data show that 
an optical group is indeed centered on NGC 6159 and contains as many as 13 members, 7 
of which are within 0.5 h^Mpc (Table 5). 
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The X-ray emission from A2199 is quite symmetric relative to other clusters (Mohr et al. 
1995, MVFS), suggesting that the inner region of A2199 has not undergone any recent major 
mergers. The bound systems are all located at projected radii significantly larger than the 
virial radius (r v « 1.6 h~ 1 Mpc, see §5.6). We identify galaxies within 1.5 /i _1 Mpc of the 
bound systems in Figure 11. NRGs385, NRGs388, and A2199 are roughly colinear, and the 
galaxy contours of A2199 are noticeably elongated along this axis (Figure 6). This alignment 
may be coincidental or it may indicate the presence of a filament of galaxies and/or dark 
matter (see Rines et al. 2001b, for further discussion). 

We next investigate the robustness of the infall mass profile with respect to variations 
in the D99 prescription. 

5.1. Choice of Smoothing Parameter 

The smoothing parameter q defines the width of the smoothing kernel in angular distance 
to the smoothing length in redshift. There is no simple a priori definition for this parameter. 
It is thus important to quantify the systematic uncertainty due to this choice (e.g., Rines et al. 
2000, 2001a). D99 showed that the choice q=25 yields accurate mass profiles for simulated 
clusters; the choices g=10 and g=50 yield similar results. RQCM find that q=5 yields a 
mass profile for the Shapley Supercluster consistent with the sum of the X-ray masses of the 
individual systems. 

Figure 12 displays the mass profile estimated from the caustics for the choices g=10, 25, 
and 50. The profiles show excellent agreement within 4 /i _1 Mpc, indicating that our results 
do not depend strongly on the choice of this parameter. Outside 4 h^Mpc, the caustics have 
negligible amplitude due to the constraint on the first derivative of the caustics (see D99 and 
§5.3). We adopt the g=25 mass profile as the standard profile to test other variations. 

5.2. Defining the Velocity Center of the Supercluster 

In the simplest case, the velocity center of the largest system in a supercluster (in 
this case, A2199) is equal to the velocity center of the entire supercluster. However, clusters 
within a particular supercluster may have substantial peculiar velocities (Bahcall et al. 1994). 

Lucey et al. (1997) estimate that the Hubble velocity with respect to the CMB at 
the distance of A2199 is cz C mb = 9190 ± 360 km s" 1 using the Fundamental Plane and 
9550 ± 400 km s -1 using D v — o. Using the same data, Hudson et al. (1997) and Gibbons 
et al. (2001) estimate 9289 ± 307 km s _1 and 9452 ± 353 km s _1 respectively. The EFAR 
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survey (Colless et al. 2001) reports czqmb = 9483±741 km s 1 using the Fundamental Plane. 
Giovanelli et al. (1998) use the Tully-Fisher relation to estimate czcmb — 9231 ±674 km s _1 . 
The IR SBF study of Jensen et al. (2001) suggests cz CM b = 9302±500 km s^ 1 for NGC 6166, 
the cD galaxy of A2199. The weighted mean yields czqmb — 9349±157 km s" 1 for the Hubble 
velocity at the distance of A2199, compared to czqmb = cz Q + 26 km s _1 = 9182 km s' 1 
from our hierarchical center. These results suggest that A2199 does not have a significant 
peculiar velocity with respect to the CMB. 

The systemic redshift cz sup of the supercluster is between 9000 and 9500 km s -1 . The 
central cluster A2199 appears to lie at a slightly lower redshift than the Great Wall. This 
offset could mean that (1) A2199 has a significant peculiar velocity with respect to the 
centroid of the supercluster potential or (2) the large-scale structure of the Great Wall 
populates only the back side of the infall region. In case (1), the redshift center of the 
supercluster potential is likely the redshift center of the Great Wall near A2199. In case (2), 
however, the supercluster potential is centered on A2199. 

Because the mass profile depends on the squared amplitude of the caustics, a small shift 
in the systemic redshift could produce a large change in the inferred mass profile. Figures 
13 and 14 show the caustics and resulting mass profiles calculated using cz sup =9000, 9156, 
and 9500 km s~ x as the redshift center of the system and smoothing parameter q = 25. The 
mass profiles for cz sup =9000 and 9156 km s _1 are in excellent agreement; the mass inferred 
for cz sup =9500 km s _1 is somewhat larger outside R3 3 /i^Mpc, but the uncertainties are 
correspondingly larger. Thus, the uncertainty in cz sup introduces < 30% uncertainty in the 
mass within 3 /i _1 Mpc; outside this radius, the systematic uncertainty introduced is at most 
a factor of 2. These uncertainties are comparable with the intrinsic systematic uncertainties 
estimated by D99 for isolated systems. 

The X-ray group NRGs385 is located within the supercluster caustics only for cz sup =9500 
km s -1 , which is the only variation which yields a significant caustic amplitude outside 
4 /i~ 1 Mpc. Indeed, if one were to sketch the locations of the caustics by eye, one would 
likely include the envelope of galaxies outside 4 /i -1 Mpc. This mass profile gives a good 
estimate of the upper bound on the total supercluster mass. At the largest radii sampled, 
the cz sup =9500 km s _1 mass profile yields a total mass of 10 15 /i _1 M Q , roughly twice the esti- 
mates of total supercluster mass within 10 h~ 1 Mpc for the standard profile. The supercluster 
mass within 8 — 10 ft, _1 Mpc is thus more uncertain than the mass within 3 — 4 h~ 1 Mpc. 
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5.3. Relaxing the Rules for Determining A(R P ) 

To reduce contamination from interloper galaxies, the method of D99 requires that 
dlnA(r) / dlnr < 1. If dlnA(r) / dlnr > 1, A(R P ) is replaced with a value such that dlnA(r) / d\wr = 
1. Qualitatively, this step prevents the caustics from "flaring out" with radius. This prop- 
erty is satisfied by all mass profiles with density profiles decreasing as p oc r~ a at large radii 
for a > 2 (all the models considered in §2.1 satisfy this requirement). This requirement 
may introduce a bias into the final mass profile. Figure 15 shows the shape of the caustics 
calculated without this requirement and without requiring that the caustics be symmetric 
(i.e., the amplitudes of the upper and lower caustics are independent). Rines et al. (2001b) 
use this variation. Figure 16 shows the mass profile calculated from these caustics using the 
minimum caustic amplitude at each radius for q = 25 and cz sup = 9156 km s -1 . Within 
4 /i~ 1 Mpc, the mass profile calculated in this manner is nearly identical to the standard pro- 
file. Outside 4 /i _1 Mpc, the mass profile lies between the standard profile and that calculated 
with cz s ,.. n =9500 km s -1 . 



5.4. Substructure and Caustics 

In A2199, we directly see the substructure which disrupts the sharp caustic pattern 
expected for simple spherical infall (Kaiser 1987). The location of subclusters within the 
caustic pattern of A2199 demonstrates that they are falling into the supercluster. We split 
the entire region into two halves along a line passing through A2199 to test the importance 
of A2197 in determining the location of the caustics. The caustics at the radius of A2197 
have a larger amplitude towards A2197 than away from it (Figure 17). There is a noticeable 
deficit of galaxies at the radial distance of A2197 on the opposite side of A2199. However, the 
caustics away from A2197 in the vicinity of NRGs388 appear similar to those towards A2197, 
although the estimator yields a smaller amplitude due to the deficit of galaxies at smaller 
radii and the constraint that the caustics not rise too quickly. We discuss the interplay of 
group and supercluster dynamics in §6. Figure 18 shows the mass profiles calculated towards 
and away A2197; the two profiles agree within uncertainties. 



5.5. Comparison to Model Mass Profiles 

We fit the mass profile to the models discussed in §2.1. Table 7 lists the best-fit param- 
eters a (fixed for SIS) and M(a) for the three models for the variations in q, cz sup , and fitting 
technique described above. The fits only include radii with A(r) > 100 km s -1 . Because 
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the data points in the mass profile are not independent, the absolute values of xl listed 
in Table 7 are only indicative, but it is clear that, as for the simulated superclusters, the 
NFW and Hernquist profiles provide acceptable fits to the infall mass profile while the SIS 
is excluded for all variations. The NFW profile generally provides a better fit to the data 
than the Hernquist profile. The parameters of the best-fit mass profile are quite robust with 
respect to the variations in the estimator discussed above. As noted in §4.2, the mass profile 
calculated assuming cz sup =9500 km s -1 yields a significantly larger total mass than the 
other estimates, but the mass profile inside 3 ft, _1 Mpc is very similar to the others. 

For the standard mass profile, the best-fit NFW parameters are log 5 C = 4.4, concen- 
tration c = r 200 /a = 8, and M 200 = 3.2 x 10 14 = 7.8M* for the NFW definition of M* 
in a ACDM model (M 200 = 20M* for the SCDM model in NFW). Our results match the 
predicted 5 C — M 20 o relation for most of the cosmological models in NFW. 

5.6. Comparison to Other Mass Estimates 

MVFS use ROSAT X-ray images and spatially-resolved spectroscopy from ASCA to 
estimate the mass profile of A2199 in the range 0.1 — 0.6 /i _1 Mpc. This mass profile is 
one of the few X-ray cluster mass profiles measured to such a large radius with a resolved 
temperature profile. The X-ray mass profile, which is completely independent of the caustic 
technique, shows excellent agreement with the standard profile (Figure 19). Figure 19 also 
shows a mass estimate from deprojecting Einstein data (White et al. 1997) and the best- 
fit NFW profile to the deprojected Chandra data at small scales (Johnstone et al. 2002, 
however, note that this NFW profile is formally rejected by the data). The standard mass 
profile agrees well with the former but not with the latter. The ICM is far from hydrostatic 
equilibrium in the core of A2199, so the mass estimate at these small radii is perhaps not 
reliable. 

We use the prescription of Girardi et al. (1998) to estimate a virial radius r v = 1.6 ft, _1 Mpc 
for A2199. We estimate the virial mass of A2199 using all galaxies inside the caustics and 
within r v . The virial mass is 5.34 x 10 14 /i~ 1 M Q , in good agreement with Girardi et al. (1998). 
The standard profile yields a mass about 30% smaller than this value at r v . This difference 
is within the uncertainty range expected for the caustic technique. 
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5.7. Velocity Dispersion Profile 

As a final consistency check of the caustic mass profile, we calculate the velocity dis- 
persion as a function of radius in bins of 25 galaxies. For this calculation, we include all 
galaxies inside the caustics of either the standard caustics or the caustics calculated with 
cz sup = 9500 km s~ x . Including only galaxies in the former caustics yields similar results 
within 4 h^Mpc; very few galaxies are classified as members beyond this radius. Figure 20 
displays both the regular and the integrated velocity dispersion profiles along with the pre- 
dicted profiles based on the best-fit Hernquist and NFW mass profiles and an assumption of 
isotropic orbits. The predicted profiles match the observed profiles quite well, and thus the 
galaxy orbits in A2199 are consistent with isotropic. Note that, as in A576, the integrated 
velocity dispersion decreases smoothly with radius and does not reach an asymptotic value 
which could be easily interpreted as the velocity dispersion of the cluster. 

6. Discussion 

Projection effects introduce the largest uncertainty in the caustic technique both for 
isolated clusters and simple superclusters (D99). The X-ray masses for A2197W and A2197E 
are factors of 6 and 13 smaller than A2199. Tests of similar systems from simulations indicate 
that the caustic technique should yield an accurate mass profile for the supercluster. The 
uncertainty in the redshift center of the supercluster results in a factor of two uncertainty in 
the total supercluster mass within 8-10 ft, _1 Mpc. Within 3 /i~ 1 Mpc, however, we estimate 
the error on the enclosed mass to be only about 30%, comparable to the intrinsic uncertainty 
estimated by simulations. 

6.1. Do Supercluster Dynamics Dominate Group Dynamics? 

At the projected radii of A2197W and A2197E, the inferred group masses are a small 
fraction of the enclosed supercluster mass. The supercluster dynamics may dominate the 
observed dynamics of groups in the infall region. The supercluster may increase the velocity 
dispersion of the groups either directly or by introducing interlopers which are bound to the 
supercluster but are not group members. 

One test of the relative importance of group and supercluster dynamics is to predict 
the velocity dispersion of the groups from their X-ray luminosities or temperatures and the 
Lx — cr or a — Tx relation (Xue & Wu 2000, MBGR and references therein). We list the 
predicted and observed velocity dispersions in Table 8. Except for A2199 and the NGC 6159 
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group, the observed velocity dispersions are larger than those predicted by the Lx — a or 
the a — Tx relation, suggesting that the observed dynamics of the groups are dominated by 
the supercluster, either physically (e.g., due to tidal forces from filaments) or observationally 
due to the increased density of interlopers. 

At the radii of A2197, NRGs388, and NGC 6159, the amplitude of the caustics is 
relatively unaffected by the presence of these subsystems. At the projected distance of 
NRGs385, however, the caustic pattern (Figure 15) and the velocity dispersion profile (Figure 
20) shows a larger amplitude than at smaller radii. This result shows that the caustic 
amplitude is strongly affected by substructure at this radius. The redshift center of the 
supercluster determines whether NRGs385 is inside the supercluster caustics. 

6.2. X-ray Contribution of Radio Sources 

X-ray emission from radio sources can contribute significantly to the total X-ray lumi- 
nosity of a galaxy group. Indeed, Feretti et al. (1995) find that both NRGs385 and NRGs388 
contain radio galaxies which are detected in X-rays. From ROSAT PSPC data, they find 
that NGC 6107 and NGC 6109 contain 6% and 1.8% respectively of the total X-ray lu- 
minosity of NRGs385. For NRGs388, they suggest that the observed X-ray luminosity is 
entirely associated with the halo of NGC 6137 and not with an intragroup medium. Reliable 
estimates of a p from the Lx — o v relation therefore require pointed X-ray observations of 
groups to determine the X-ray luminosity due to the intragroup medium. 

6.3. Are the Groups Really Infalling? 

The location of the groups within the caustics indicates that they are dynamically 
linked to the A2199 supercluster. This connection, however, does not necessarily imply that 
the groups are currently infalling onto the supercluster. The turnaround radius r ta for a 
A = universe satisfies p m (< r ta ) = Q m p c (l + A 0tturn ) where p c is the critical density and 
Aojum is given by Equation 8 in Regos & Geller (1989). In the limit Q m = 1, p m (< r ta) = 
(97r 2 /16)p c . The standard mass profile yields r ta ?^4.5-5.7 /i _1 Mpc for Q m — 1 — > 0.025 
where we take a lower limit on Q m from the baryon density Qbh 2 = 0.025 (Pettini & Bowen 
2001) and the assumption that h < 1. For the more massive profile inferred by assuming 
cz sup =9500 km s -1 , the turnaround radius for the above range of Vt m is r ta =5.0-6.5 /i _1 Mpc. 
The turnaround radius is smaller in the presence of a non-zero cosmo logical constant, but 
this effect is small in the present epoch for A < 1 (Lilje & Lahav 1991; Martel 1991). Thus, 
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A2197W, A2197E, and NRGs388 are in the infall region of A2199 for most combinations of 
geometry and cosmology (the expectation value of the deprojected radius is < r >= nR p /2). 
The NGC 6159 group is a borderline case; it is infalling for the standard mass profile and 
our adopted cosmology only for small angles (< 20°) between the plane of the sky and the 
line connecting A2199 and NGC 6159. NRGs385 and NRGs396 are not infalling onto the 
supercluster. A simulation of the future evolution of the nearby universe (Nagamine & Loeb 
2002) shows that some objects outside the turnaround radius of clusters in the present epoch 
but with large peculiar velocities may eventually become cluster members. 

Based on the presence of a caustic pattern around NRGs385 and an unusual offset be- 
tween its X-ray and optical peaks, we previously suggested that NRGs385 might be infalling 
onto A2199 (Rines et al. 2001b). We have shown here, however, that the mass profile ro- 
bustly excludes this possibility. The unusual morphology of NRGs385 is probably related to 
the local large-scale structure and is not caused by effects related to infall (e.g., ram pressure 
stripping). 

The caustic pattern is present beyond the turnaround radius, suggesting that the gravi- 
tational potential of a cluster (or simple supercluster) and surrounding large-scale structure 
can significantly affect the kinematics of objects outside the infall region. The caustics ex- 
tend beyond the turnaround radius in the simulations of D99. A similar effect is seen in the 
galaxy-galaxy autocorrelation function (e.g., Peacock et al. 2001). The existence of large- 
scale structure obscures the location of the turnaround radius in the projected radius-redshift 
diagram. 

7. Conclusions 

We have previously demonstrated that the infall region of A2199 contains several bound 
subclusters (Rines et al. 2001b). Here, we show that, in simulations, caustics yield accurate 
mass profiles of superclusters dominated by a central cluster 2-6 times more massive than 
all other systems in the supercluster. The mass profiles of the simulated clusters agree with 
NFW or Hernquist models but not with a singular isothermal sphere. 

X-ray data demonstrate that A2197 is the superposition of two bound X-ray groups (cen- 
tered on the bright galaxies NGC 6160 and NGC 6173) which may be interacting. Because of 
the complexity of A2197, the virial theorem does not yield an accurate mass estimate. The 
X-ray masses for A2197W and A2197E are factors of 6 and 13 smaller than A2199. Tests 
of similar superclusters from simulations indicate that the caustic technique should yield an 
accurate mass profile for the A2199 supercluster. 
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The supercluster mass profile suffers some uncertainty due to the complexity of the 
surrounding large-scale structure and the resulting difficulty in determining the redshift 
center of the supercluster. Within 3 /i -1 Mpc, however, this uncertainty is no larger than 
the scatter predicted by numerical simulations for isolated clusters. The infall mass profile 
is in excellent agreement with independent X-ray estimates at small radii. At large radii, 
the mass profile is clearly inconsistent with an isothermal sphere but agrees with profiles 
motivated by simulations. 

The mass profile indicates that A2197W, A2197E, and NRGs388 are bound, infalling 
groups; NRGs385 and NRGs396 (and the NGC 6159 group for inclination angles > 20°) are 
not infalling onto the supercluster in the present epoch. This result demonstrates that the 
presence of caustics at large radii from clusters does not prove that galaxies or groups within 
these caustics are inside the infall region. In other words, the gravitational potential of a 
cluster and surrounding large-scale structure can significantly affect the kinematics of objects 
outside the turnaround radius. The observed velocity dispersions of the groups surrounding 
A2199 are significantly larger than predicted by their X-ray properties and scaling relations. 
The observed kinematic properties of galaxy groups within and near superclusters are likely 
dominated by the supercluster and not by the group potential. 

The caustic technique is a robust estimator of cluster mass profiles ranging in mass from 
Fornax to the Shapley Supercluster. We demonstrate that the technique is robust for simple 
superclusters, both in simulations and in observations of the A2199/A2197 supercluster. We 
plan to analyze a sample of 8 nearby clusters and study the mass-to-light ratio throughout 
the infall regions of several of these systems. 
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Fig. 1. — Redshift versus radius and the caustics for a simulated supercluster (ID 121) with 
mass ratio Mi/M 2 = 2.0 (see text for definition) projected along three different lines of sight. 
We use different line styles to enable comparison with the estimated mass profiles in Figures 
4 and 5. 
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Fig. 2. — Same as Figure 1 for the simulated supercluster ID 156 with mass ratio Mi/M 2 = 
AA. 
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Fig. 3. — Same as Figure 1 for the simulated supercluster ID 162 with mass ratio Mi 
6.2. 
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Fig. 4. — Comparison of caustic mass profiles to true mass profiles for simulated superclus- 
ters. The three sets of lines show the 1-a ranges of the infall mass profile from three different 
lines of sight. The line styles are the same as in Figures 1-3. 
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Fig. 5. — Difference between estimated and true masses of the simulated superclusters in 
units of the uncertainty in the estimated mass. The estimates differ by more than 2a from 
the true mass profile only in the supercluster with Mi/M 2 = 2.0. The line styles are the 
same as in Figures 1-3. 
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Fig. 6. — Galaxy density contours overlaid on RASS X-ray data. The RASS data have been 
exposure corrected. Both RASS data and the galaxy density were smoothed with a Gaussian 
with 2' FWHM. Some groups of galaxies have no associated X-ray emission detectable in the 
RASS. The inner circle is the limit of our redshift survey; our sample also includes galaxies 
with redshifts from the literature within the outer circle. 




Fig. 7. — X-ray contours from a 13.5-ksec ROSAT PSPC observation of A2197 overlaid on 
optical POSS I data. The diameter of the PSPC field is 2°. North is up and East is to the 
left. 
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Fig. 8. — Redshift histogram of galaxies within 1.5 h x Mpc of A2197. The solid line shows 
the best-fit Gaussian and the dashed line shows a two-component model. Arrows indicate the 
redshifts of NGC 6173 and NGC 6160, the cD galaxies of A2197E and A2197W respectively. 
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Fig. 9.— ASCA X-ray spectrum of A2197W. The data points show the GIS 2 and GIS 3 
instruments and the lines show the best-fit thermal plasma model convolved with the ASCA 
response. 
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Fig. 10.— Same as Figure 9 but for A2197E. 
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Fig. 11. — Projected radius versus redshift for galaxies surrounding A2199. Lines indicate 
our estimate of the caustics. Crosses, hexagons, open squares, asterisks, and pentagons 
indicate galaxies in A2197, NRGs385, NRGs388, NRGs396, and NGC6159 respectively. 
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Fig. 12. — Mass profile of the A2199 supercluster. The thick solid line indicates the mass 
profile estimated using q = 25 (see text for definition), the thin solid lines and shaded region 
show the l-o" uncertainties in this profile. The sets of dashed and dash-dotted lines indicate 
the mass profiles and 1-a uncertainties for q=10 and 50 respectively. 
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Fig. 13. — Same as Figure 11, but varying the systemic redshift of the supercluster cz sup . 
Lines indicate our estimate of the caustics. Dashed, solid, and dash-dotted lines show the 
caustics calculated with cz sup = 9000, 9156, and 9500 km s _1 respectively. 
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Fig. 14. — Effect on the mass profile of varying the central redshift of the supercluster. 
Dashed, solid, and dash-dotted lines show the caustics calculated with cz sup = 9000, 9156, 
and 9500 km s _1 respectively. 
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Fig. 15. — Same as Figure 11, but using relaxed requirements to calculate the caustics. 
Dashed lines show the caustics on either side of the cluster; the solid lines show the minimum 
caustic amplitude at each radius. The latter caustics are defined primarily by the lower 
caustic. 
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Fig. 16. — Mass profile calculated with relaxed requirements on the slope of the caustics. 
The dashed lines show the standard caustics for reference. 
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Fig. 17. — Projected radius versus redshift for galaxies surrounding A2199. The top panel 
shows galaxies towards A2197; the bottom panel shows galaxies away from A2197. Dashed 
lines indicate the caustics calculated for each sample; the solid lines are the caustics from 
the full sample. Crosses indicate galaxies in A2197, NGC 6159, and NRGs396 (top panel) 
and NRGs385 and NRGs388 (bottom panel). 
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Fig. 18. — Mass profile calculated towards (solid lines) and away (dashed lines) from A2197. 
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Fig. 19. — Comparison of infall mass profile to independent X-ray estimates. The 1-er range 
of the infall mass profile is shown in red and the 90% confidence range of the X-ray mass 
profile is shown in blue. The NFW profile which best fits the X-ray data is indicated by a 
thick solid line. The point indicates a deprojected X-ray mass estimate and the dashed line 
shows the NFW profile which best fits Chandra X-ray data. 
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Fig. 20. — Velocity dispersion profile of A2f99 member galaxies (open squares). Filled 
squares show the integrated velocity dispersion profile. Dotted and dot-dashed lines show 
the profiles predicted by the Hernquist and NFW profiles respectively assuming isotropic 
orbits. The predicted profiles have no free parameters. 
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Table 1: Properti es of Simulated Simple Superclusters 

System Mi Mi/M 2 M subhalo R su bhaio 

1O 13 /i" 1 M lO^h^M® h^Mpc 



121 


18.2 


1.95 


9.3 


3.19 








8.8 


3.55 








4.1 


4.96 








4.7 


5.87 


156 


33.5 


4.43 


7.6 


2.90 








2.1 


3.66 








1.2 


6.74 








2.5 


10.13 


162 


31.0 


6.20 


3.2 


2.74 








5.0 


3.38 








1.1 


7.97 








2.0 


8.89 



Table 2: Simulated Mass Profile Fit Parameters 



System 


Profile 


a 


M(a) 
1O 14 M 


X 2 


V 


121 


NFW 


0.24±0.01 


0.38±0.02 


101 


23 


121 


Hcrnquist 


0.49±0.02 


0.98±0.05 


167 


23 


121 


SIS 


0.5 


0.42±0.01 


1152 


24 


156 


NFW 


0.18±0.01 


0.48±0.02 


18 


23 


156 


Hernquist 


0.41±0.02 


1.44±0.07 


64 


23 


156 


SIS 


0.5 


0.93±0.03 


790 


24 


162 


NFW 


0.25±0.01 


0.71±0.04 


52 


23 


162 


Hernquist 


0.53±0.02 


1.93±0.10 


78 


23 


162 


SIS 


0.5 


0.69±0.02 


1128 


24 
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Table 3. Spectroscopic Data for A2199 a 



RA 


DEC 


cz 




Reference 


(J2000) 


(J2000) 


(km s _1 ) 


(km s _1 ) 




15 57 54.36 


+41 56 13.8 


10387 


125 


04 


15 57 58.01 


+41 47 32.6 


10464 


033 


04 


15 57 59.40 


+40 02 02.0 


09300 


000 


04 


15 58 47.35 


+41 56 17.2 


11561 


033 


04 


15 59 28.85 


+39 49 39.5 


09296 


045 


04 



a The complete version of this table is in the electronic edition 
of the Journal. The printed edition contains only a sample. 

References. (1) FAST spectra (2) Mohr et al. 2002, in 
prep.; (3) Mahdavi et al. 1999; (4) Falco et al. 1999; (5) Kirsh- 
ner et al. 1983: (6) Strauss & Huchra 1988; (7) de Vaucouleurs 
et al. 1991; (8) Freudling et al. 1992; (9) Zabludoff et al. 1993; 
(10) Haynes et al. 1997; (11) Hill & Oegerle 1998. 
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Table 4: ARCHIVAL POINTED X-RAY OBSERVATIONS 



System 


Satellite 


Instrument 


Sequence 


Live Time(ks) 


A2199 


ASCA 


GIS/SIS 


80023000 


67 




ROSAT 


PSPC 


rp800644n00 


41 


A2197W 


ROSAT 


PSPC 


rp800363n00/a01 


13.5 




ASCA 


GIS/SIS 


85068000 


45 


A2197E 


ROSAT 


PSPC 


rp800363n00/a01 


13.5 




ASCA 


GIS/SIS 


85069000 


49 


NRGs385 


ROSAT 


PSPC 


rp700508n00/a01 


5.6 



Table 5: PROPERTIES OF SYSTEMS ASSOCIATED WITH A2199 



System 


X-ray Coordinates 


cz 




R p 


logLx 


N 3 * p 


N tot 




RA (J2000) 


DEC (J2000) 


km s _1 


km s _1 


ft _1 Mpc 


h~ 2 ergs s _1 






A2199 


16 28 38 


39 33 05 


9101±50 


796+ 33 




44.1 


255 


339 


A2197W 


16 27 41 


40 55 40 


9144±52 


584+^0 


2.2 


42.5 


128 


144 


A2197E 


16 29 43 


40 49 12 


9100±53 


595+g 


2.0 


42.4 


126 


151 


NRGs385 


16 17 15 


34 55 00 


9308±86 


643+^ 


8.1 


42.8 


59 


66 


NRGs388 


16 23 01 


37 55 21 


9421±74 


5631S 


3.1 


42.3 


58 


78 


NRGs396 


16 36 50 


44 13 00 


9554±128 


513173 8 


7.7 


a 

42.0 


16 


29 


NGC 6159 


16 27 25 


42 40 27 


9566±95 


344l^°° 


4.9 


b 

42.3 


13 


22 



"Detection significance = 2.7cr 

fc HRI observation (Trinchieri & Pietsch 2000) 
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Table 6: Properties of X-ray Sources in A2197 field with Optical Counter- 
parts 



ID 


Type 


X-ray Coordinates (J2000) 


cz 


fx (0. 1-2.4 keV) 


\ogL x 






RA 


DEC 


km s _1 


10~ 12 ergs cm~ 2 s _1 


h~ 2 ergs s _1 


A2197W 


Group 


16 27 41 


40 55 40 


9408±77 a 


3.13±0.10 


42.5 


A2197E 


Group 


16 29 43 


40 49 12 


8842±45 a 


2.57±0.09 


42.4 


NGC 6146 


Galaxy 


16 25 10 


40 53 34 


8879±25 


0.18±0.03 


41.3 


NGC 6150 b 


Galaxy 


16 25 50 


40 29 19 


8707±31 


0.04±0.01 


40.6 


NGC 6175 


Galaxy 


16 29 58 


40 37 43 


8986±34 


0.60±0.04 


41.8 


III Zw 078 b 


Galaxy 


16 26 46 


40 41 40 


8838±40 


0.07±0.01 


40.9 


Anonymous 


Galaxy 


16 28 26 


40 51 39 


9976±24 


0.10±0.01 


41.0 


2MASXi J1626069+412046 


Galaxy 


16 26 07 


41 20 39 


17592±21 


0.05±0.01 


41.3 


A2187 


Cluster 


16 24 11 


41 13 46 


54938±54 a 


2.37±0.08 


44.2 


A2196 


Cluster 


16 27 26 


41 29 55 


38811±30 a 


1.28±0.06 


43.6 


FBQS J162901.3+400759 


QSO 


16 29 01 


40 08 34 


81545 d 


1.41±0.06 


44.3 


KUV 16295+4054 


QSO 


16 31 13 


40 48 40 


77048 6 


1.54±0.08 


44.5 



a Redshift of bright galaxy near center 

^tentative identification, X-ray-optical separation ~ 1' 

c see also Knezek & Bregman (1998) 

d Bade et al. (1995) 

e Crampton et al. (1992) 
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Table 7: Mass Profile Fit Parameters 



1V1L, l ..in. J 1.1 


T-*r , f~\fi 1 r> 

JT I UIIIL- 


Q 


CI 


1 — rr 




1 — rr 


V 2 

A 










/i -1 Mpc 




10 14 M Q 


1O 14 M 






^ T o i~i /~\ ovrl 
O LdllCld-I (_1 


in r vv 




1 7 






n A^ ^8 


1 9 

J. . Zi 


9zL 

Z4: 


cm H tiyH 
O LdjIICIdl LI 


1— 1 nrn nmcf 
11L.I IIL[U.Ici I 


1 n 


U. OK) 






1 fin i on 


1 1 

j. . j. 




Sl+ 'A nrl a rH 

kJ LOiIlUQil 


SIS 


10 


0.5 




1.01 


97-1 06 


36.1 


25 


OLdllQcll Q 


in r vv 


9^ 


1 A 


nil n 1 q 


A7 


n ad n V\ 

u.^tu-u.oo 


n ^ 


9^, 


DlcLIlUdrl Q 


n.L.1 I1L[ Ulo L 


9^, 


A7 


u.oy-u.o / 


l.UO 


1 m 1 Qf; 

1 .OU-l.Od 


1 9 


9^, 

zo 


Strinrlrirr] 


SIS 


25 


0.5 




1.00 


94-1 05 


28.8 


26 


Standard 


NFW 


50 


0.15 


0.10-0.23 


0.49 


0.36-0.65 


0.2 


23 


Standard 


Hernquist 


50 


0.46 


0.33-0.64 


1.65 


1.35-2.05 


0.4 


23 


Standard 


SIS 


50 


0.5 




1.08 


0.96-1.18 


7.5 


24 


cz sup = 9000 


NFW 


25 


0.13 


0.10-0.15 


0.43 


0.37-0.48 


1.7 


23 


cz sup = 9000 


Hernquist 


25 


0.42 


0.36-0.49 


1.53 


1.40-1.66 


1.5 


23 


cz sup = 9000 


SIS 


25 


0.5 




0.97 


0.93-1.02 


52.0 


24 


cz sup = 9500 


NFW 


25 


0.30 


0.26-0.34 


0.79 


0.73-0.86 


2.1 


61 


c^ s «p = 9500 


Hernquist 


25 


0.94 


0.86-1.04 


2.80 


2.62-2.96 


10.2 


61 


cz sup = 9500 


SIS 


25 


0.5 




0.81 


0.78-0.85 


121 


62 


Relaxed 


NFW 


25 


0.15 


0.11-0.19 


0.46 


0.41-0.52 


0.6 


54 


Relaxed 


Hernquist 


25 


0.63 


0.52-0.75 


1.96 


1.82-2.11 


2.5 


54 


Relaxed 


SIS 


25 


0.5 




0.59 


0.56-0.62 


69.0 


55 



Table 8: PREDICTED AND OBSERVED VELOCITY DISPERSIONS 



System 


L x 


Tx 




a P,T x 


Up 


A2199 


44.1 


4.5 


892 


810 


796+jg 


A2197W 


42.5 


1.55 


357 


423 


584 +4 o 


A2197E 


42.4 


0.96 


337 


316 


5951S 


NRGs385 


42.8 




424 




643+^ 


NRGs388 


42.3 




318 




563tg 


NRGs396 


42.0 




268 




513 + 7f 


NGC6159 


42.3 




318 




344+™° 



